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Electrokinetic  power  generation  efficiency  using  a  two-dimensional  axisymmetrical  model  is  numerically 
investigated.  A  finite-length  nanoscale  surface-charged  cylindrical  capillary  with  reservoirs  connected 
at  the  capillary  ends  is  considered  as  the  physical  domain.  The  Navier-Stokes,  Laplace,  Poisson,  and 
Nernst-Planck  equations  are  solved  simultaneously  to  obtain  the  fluid  flow,  electrical  potential,  ion  con¬ 
centration  and  electrical  current  in  the  flow  field.  The  energy  conversion  efficiency  predicted  using  a 
one-dimensional  model  assuming  an  infinitely  long  channel,  Boltzmann  ion  distribution  and  equal  ionic 
electrical  mobility  is  also  carried  out  and  compared  with  the  two-dimensional  result. 

The  two-dimensional  model  results  show  that  the  electrostatic  potential  gradient  resulting  from  the 
concentration  changes  at  the  capillary  entrance  and  exit  and  fluid  flow  produce  a  conductive  current 
that  reduces  the  total  current  in  the  flow  field.  The  conductive  current  due  to  the  electrostatic  potential 
gradient  increases  with  the  decrease  in  electrolyte  bulk  concentration  and  increase  in  surface  charge 
density.  This  results  in  nonlinear  variations  in  the  electric  current-voltage  curve  and  maximum  con¬ 
version  efficiency  as  functions  of  the  surface  charge  density  and  dimensionless  Debye  length  when  the 
electrolyte  bulk  concentration  is  low.  Comparison  of  the  maximum  efficiencies  predicted  from  one-  and 
two-dimensional  models  indicates  that  the  one-dimensional  model  is  valid  only  when  the  dimensionless 
Debye  length  is  large  and  the  surface  charge  density  is  small  because  the  electrostatic  potential  gradient 
is  neglected.  The  two-dimensional  model  also  predicts  that  optimum  maximum  energy  conversion  effi¬ 
ciency  can  be  obtained  when  the  dimensionless  Debye  length  is  equal  to  2  and  its  magnitude  increases 
with  the  increase  in  surface  charge  density. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Because  of  the  developments  in  recent  microfabrication  tech¬ 
nology,  micro  or  nanoscale  channels  have  become  available,  playing 
an  important  role  in  many  engineering  applications.  For  example, 
fluid  transport  in  nanoscale  fluidic  systems  by  applying  an  exter¬ 
nal  electrical  field  is  used  for  chemical  and  biological  analyses  in 
the  lab-on-a-chip  devices  [1-3].  Besides  chemical  and  biological 
analyses,  micro  and  nanoscale  fluidic  system  applications  in  energy 
conversion  based  on  electrokinetic  effects  have  recently  gained 
much  interest  in  the  research  society  because  of  the  demand  for 
finding  new  and  environmental-friendly  energy  resources  [4,5]. 

The  electrokinetic  energy  conversion  principle  utilizes  two 
fundamental  electrokinetic  effects:  the  streaming  potential  and 
electroosmosis.  These  two  phenomena  rely  on  the  existence  of 
an  electrical  double  layer  (EDL)  near  a  channel  wall.  The  theo¬ 
ries  regarding  the  formation  of  EDL,  electroosmosis  and  streaming 
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potential  are  well  documented  in  several  textbooks  [6-8].  This  is 
based  on  the  fact  that  most  surfaces  are  electrically  charged  when 
they  are  brought  into  contact  with  a  polar  medium  such  as  an 
aqueous  electrolyte.  Because  of  the  charged  surface,  preferentially 
distributed  counter-ions  and  co-ions  are  created  near  the  surface. 
The  ion  distribution  combined  with  random  thermal  motion  results 
in  EDL  formation.  As  the  fluid  flow  is  driven  by  an  externally  applied 
pressure  gradient,  the  fluid  flow  carries  the  ions  to  the  channel 
downstream.  The  charge  transport  results  in  an  electric  current, 
known  as  the  streaming  current.  The  charge  polarization  due  to 
accumulated  ions  at  the  channel  downstream  creates  an  electrical 
potential  difference  between  the  channel  ends  and  a  conductive 
current  with  a  direction  opposite  to  the  flow  direction.  The  electri¬ 
cal  potential  difference  is  termed  the  streaming  potential  because 
it  is  created  by  the  fluid  flow.  The  streaming  potential  theory  has 
recently  been  applied  to  build  an  electric  power  generation  unit 
since  the  maximum  streaming  potential  can  be  regarded  as  an 
open  circuit  potential  (or  the  electromotive  force)  similar  to  that 
in  an  electrochemical  cell  [9-14].  Because  of  this  similarity,  electric 
power  generation  based  on  the  streaming  potential  theory  may  be 
termed  the  “Electrokinetic  Cell”  (EI<  cell). 
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Nomenclature 

a  radius  of  cylindrical  capillary  (m) 

AR  radius  of  computational  domain  in  reservoirs  (m) 
c0  bulk  electrolyte  molar  concentration,  M  (mol  L-1 ) 

D  molecular  diffusivity  (m2  s-1 ) 

F  Faraday  constant  (96,500  C  mol-1) 

G  hydrodynamic  conductance  (m3  Pa-1  s-1 ) 

L  length  of  cylindrical  capillary  (m) 

Lr  length  of  computation  domains  in  reservoirs  (m) 

m  electrical  mobility  (mol  s  kg-1 ) 

M  streaming  conductance  (S) 

N  number  of  ions  in  the  flow  field 

p  pressure  (Pa) 

Pin  input  pumping  power  (W) 

P out  output  electric  power  (W) 

Q  flow  rate  (m3  s-1 ) 

R  universal  gas  constant  (8.314J  mol-1 1<-1 ) 

5  electrical  conductivity  (S) 

T  fluid  temperature  (I<) 

V  fluid  flow  velocity  vector 

z  valence  of  ion 

Z  figure  of  merit  of  energy  conversion 

Greek  letters 

A  electrolyte  electrical  conductivity  (S  m-1 ) 

s  dielectric  permittivity, 

£ o  permittivity  of  vacuum  (8.854  x  10-12  C  V-1  m-1 ) 

0  externally  applied  voltage  (V) 

r\  energy  conversion  efficiency 

fi  fluid  viscosity  (Pa  s) 

pe  net  charge  density  (Cm-3 ) 

a  surface  charge  density  (C  m-2 ) 

0  electrostatic  potential  due  to  surface  charge  (V) 

Subscript 

i  ith  ion 

max  maximum  value 

1  inlet  reservoir 

2  outlet  reservoir 


From  thermodynamic  point  of  view,  the  EK  cell  can  be  viewed 
as  an  energy  conversion  process.  The  input  power  is  the  pumping 
power  (product  of  pressure  difference  applied  and  flow  rate)  and 
the  output  is  the  electrical  power  (product  of  streaming  potential 
created  and  electric  current).  A  pioneering  study  dealing  with  elec¬ 
tro  kinetic  energy  conversion  is  the  work  of  Osterle  [15].  Based  on 
Osterle’s  work  one-dimensional  analysis  of  electrokinetic  energy 
conversion  thermodynamic  efficiency  has  been  reported  in  sev¬ 
eral  studies  [13-14,16].  The  one-dimensional  model  is  based  on  the 
infinitely  long  channel  length,  equal  ion  diffusivity  and  mobility, 
and  Boltzmann  ion  distribution  assumptions.  These  assumptions 
may  lead  to  under  or  over-predicted  energy  conversion  efficiency. 

In  order  to  raise  the  energy  conversion  efficiency  using  the  elec¬ 
trokinetic  effect,  a  flow  field  containing  many  flow  passages  must  be 
employed  to  increase  the  electric  current.  These  multiple  channel 
structures  can  be  obtained  through  microfabrication  technology, 
chemical  synthesis  or  natural  structures  [9-10,17].  In  general,  these 
microstructures  can  be  viewed  as  a  charged  porous  media  as  they 
are  in  contact  with  a  polar  medium.  Because  of  the  structural 
characteristics,  the  charged  porous  media  can  be  considered  as 
composed  of  an  array  of  identical  parallel  charged  cylindrical  capil¬ 
laries  with  micro  or  nanoscale  size.  The  charge  inside  the  porous 
media  can  be  considered  uniformly  distributed  on  the  capillary 


walls  with  a  certain  charge  density.  From  this  point  of  view,  it  is 
more  appropriate  to  consider  the  channel  in  the  energy  conver¬ 
sion  unit  as  a  charged  cylindrical  capillary  instead  of  a  planar  slit  or 
rectangular  channel,  as  investigated  in  the  studies  of  Daiguji  et  al. 
[11-12]  and  van  der  Fleyden  et  al.  [13-14]. 

In  this  study  a  two-dimensional  EK  cell  model  is  established 
and  solved  numerically.  A  nanoscale  finite-length  surface-charged 
cylindrical  capillary  including  the  inlet  and  outlet  reservoirs  is 
considered  as  the  physical  domain  [11-12,18].  The  fundamental 
equations  governing  the  fluid  flow,  multi-ion  transport  and  elec¬ 
trical  potential  distribution  are  solved  simultaneously  without 
the  assumptions  made  in  the  one-dimensional  model.  The  ther¬ 
modynamic  efficiency  results  obtained  from  both  the  one-  and 
two-dimensional  models  are  compared  and  discussed. 

2.  Physical  and  mathematical  models 

As  shown  in  Fig.  1 ,  we  consider  a  cylindrical  capillary  having  a 
radius  of  a  and  a  length  of  L.  Without  loss  of  generality,  the  capillary 
wall  is  considered  negatively  charged  with  a  surface  charge  density 
of  a.  The  ends  of  the  capillary  are  connected  to  inlet  and  outlet 
reservoirs  which  are  assumed  to  be  infinitely  large  compared  with 
the  capillary  size.  We  assume  that  an  aqueous  salt  solution  with 
bulk  molar  concentration  c0  entirely  fills  the  system  and  the  system 
is  initially  in  equilibrium.  The  physical  domain  considered  in  this 
study  is  the  same  as  that  in  the  studies  of  Daiguji  et  al.  [11-12] 
except  that  the  planar  slit  is  replaced  by  a  cylindrical  capillary. 

For  the  EK  cell  operation,  fluid  flow  is  driven  by  an  applied 
pressure  difference  and  creates  an  electrical  potential  difference. 
This  can  be  described  using  the  pressure  and  voltage  differences 
between  the  inlet  and  outlet  reservoirs.  As  shown  in  Fig.  1,  pres¬ 
sure  and  potential  in  the  inlet  and  outlet  reservoirs  are  denoted 
as  0i,  pi,  02,  P2>  respectively.  The  externally  applied  pressure  dif¬ 
ference  is  Ap=p2  -p\  <0  and  the  created  streaming  potential  is 
A0  =  0 2  -0i  >0.  Under  such  applied  pressure  difference,  the  fluid 
flow  direction  is  from  the  inlet  reservoir  to  outlet  reservoir. 

2.1.  One-dimensional  model  ofEK  cell 

The  one-dimensional  EK  cell  model  has  been  given  in  several 
investigations  [13-14,16].  A  brief  summary  based  on  the  results 
reported  by  Xuan  and  Li  [16]  is  given  below.  Based  on  the  Onsager 
reciprocal  theorem,  the  relations  between  the  flow  rate  Q,  electric 
current  /,  pressure  difference  A p,  and  the  voltage  difference  A 0 
for  the  flow  with  electrokinetic  effect  through  a  channel  with  any 
cross-sectional  shape  can  be  written  as, 

Q  =  G(-Ap)  +  M(- A0)  (1) 

/  =  M(- Ap)  +  S(  — A0)  (2) 

where  G  represents  the  hydrodynamic  conductance,  M  charac¬ 
terizes  the  electroosmotic  flow  rate  in  Eq.  (1)  and  the  streaming 
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Fig.  1.  Schematic  of  the  physical  system. 
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current  in  Eq.  (2),  and  S  represents  the  electrical  conductance.  From 
the  energy  conversion  point  of  view,  the  EK  cell  can  be  treated  as 
a  device  that  converts  the  pumping  power  into  electrical  power. 
Referring  to  Fig.  1,  the  input  pumping  power,  generated  electrical 
power  and  energy  conversion  efficiency  can  be  written  as, 

Pin  =  Q(-Ap),  P0ut  =  /(A</>),  (3) 

1  in 

As  derived  in  the  study  of  Xuan  and  Li  [16],  the  maximum  efficiency 
of  EK  cell  is  given  as, 


1)  max  — 


(l  -TT^Z)2 

z 


(4) 


where  Z  is  called  the  figure  of  merit  of  electrokinetic  energy  con¬ 
version  defined  as, 


As  indicated  in  the  study  of  Xuan  and  Li  [16],  Z  depends  on  the 
dimensionless  Debye  length,  zeta  potential  and  a  parameter  called 
the  Levine  number  which  indicates  the  nominal  ratio  of  convective 
to  conductive  electrical  current  [19].  It  can  be  said,  based  on  Eqs. 
(4)  and  (5),  that  the  maximum  EK  cell  efficiency  depends  on  the 
electrochemical  properties  of  the  working  fluid  and  wall  material 
such  as  the  bulk  concentration,  electrical  conductivity,  dielectric 
constant  and  zeta  potential  for  flow  in  a  given  channel.  The  maxi¬ 
mum  efficiency  can  be  obtained  for  any  applied  pressure  difference 
since  it  is  independent  of  the  pressure  and  voltage  differences.  As 
noted  earlier,  Eqs.  (1),  (2)  and  (4)  are  derived  under  assumptions 
including  an  infinitely  long  channel,  equal  ionic  electrical  mobility 
and  Boltzmann  ion  distribution.  These  assumptions  are  removed  in 
the  following  considering  a  two-dimensional  EK  cell  model. 


2.2.  Two-dimensional  EI<  cell  model 


By  assuming  that  the  steady  state  flow  and  continuum  assump¬ 
tion  are  valid,  the  governing  equation  describing  the  fluid  flow  is 
the  modified  Navier-Stokes  equation  written  as, 

V  •  V  =  0  (6a) 

V  W  =  -Vp  +  piV2V  -  peV{<p  +  \jf)  (6b) 

where  V  is  the  fluid  flow  velocity  vector,  p  is  the  pressure  and  /n 
is  the  viscosity.  The  last  term  on  the  right-hand  side  of  Eq.  (6b)  is 
known  as  the  electric  body  force  per  unit  volume  when  electro- 
striction  and  the  permittivity  gradient  in  the  electrical  field  are 
ignored.  pe  is  the  net  charge  density  defined  as, 


where  D2  and  m2  are  the  molecular  diffusivity  and  electrical  mobility 
of  the  ith  ion  presented  in  the  system,  respectively.  For  a  dilute  elec¬ 
trolyte,  the  ion  electrical  mobility  can  be  related  to  the  molecular 
diffusivity  by  the  Nernst-Einstein  equation, 

m'=S  (io) 

where  R  is  the  universal  gas  constant  and  T  is  the  fluid  tempera¬ 
ture.  Noting  that  Eqs.  (6)-(9)  are  the  general  governing  equations 
describing  the  fluid  flow,  ion  concentration  and  potential  distribu¬ 
tions  without  any  assumptions  such  as  Boltzmann  distribution  and 
equal  ion  mobility. 

The  current  flow  in  the  flow  field  must  be  computed  to  evaluate 
the  EK  cell  efficiency.  By  knowing  the  fluid  velocity,  ion  concentra¬ 
tion  and  electrical  potential  distributions,  the  current  density  in  the 
flow  field  can  be  calculated  as, 

?  =  -  A  V(</>  +  VO  —  F  X )  ZiDi  VCi  +  ZiCi  (1! ) 

where  A  =  F2  ^zfmiq  is  the  electrolyte  electrical  conductivity.  By 
integrating  Eq.  (11)  through  the  capillary  cross-section,  the  total 
current  at  any  cross-section  can  be  obtained.  In  Eq.  (11),  the  terms 
on  the  right-hand  side  are  conductive,  diffusive,  and  the  streaming 
currents,  respectively. 

All  governing  equations  were  written  in  a  cylindrical  coordi¬ 
nate  system  (r,0,z).  Utilizing  the  symmetry  in  the  0  coordinate,  we 
can  recast  the  problem  as  an  axisymmetric  two-dimensional  model 
(r,z).  The  fluid  velocity  components  in  the  r-  and  z-directions  are  u 
and  v,  respectively.  The  boundary  conditions  are  required  to  solve 
Eqs.  (6)-(9).  Because  of  the  geometrical  symmetry,  the  computa¬ 
tional  domain  is  taken  as  half  of  the  physical  domain,  as  shown  in 
Fig.  2.  Since  the  inlet  and  outlet  reservoirs  are  relatively  large  com¬ 
pared  with  the  capillary  size,  only  portions  of  the  reservoirs  near  the 
capillary  entrance  and  exit  are  included  in  the  computation  domain. 
As  indicated  in  Fig.  2,  the  radius  and  length  of  the  computation 
domains  in  the  reservoirs  are  taken  to  be  AR  and  Lr,  respectively.  As 
indicated  in  the  studies  of  Daiguji  et  al.  [11-12]  and  Mansouri  et  al. 
[18],  as  long  as  the  computational  domain  size  in  the  reservoirs  is 
large  enough,  no  significant  effects  in  the  fluid  flow,  ion  transport 
and  potential  distribution  will  result.  Referring  to  Fig.  2  and  taking 
the  outlet  reservoir  as  the  pressure  and  potential  references,  the 
boundary  conditions  for  simulating  the  fluid  flow  velocity,  ion  con¬ 
centration  and  potential  distributions  in  the  EK  cell  are  specified  as 
follows: 

Boundary  AB: 


N 

Pe  =  F^ZjC,- 
1  =  1 


Pi  =  -Ap,  \f/  =  0,  (p-[  =  -Acj),  q  =  c0  (12a) 

Boundaries  BC,  FG  and  AH: 


where  F  is  the  Faraday  constant,  N  is  the  number  of  ions  in  the  flow 
field,  Cj  and  z2  are  the  molar  concentration  and  the  valence  of  the  ith 
ion,  respectively.  The  potential  gradient  v0  results  from  the  poten¬ 
tial  difference  between  the  inlet  and  outlet  reservoirs,  while  the 
potential  gradient  results  from  the  electrostatic  potential  dis¬ 
tribution  due  to  the  net  charge  density  in  the  flow  field.  These  two 
potential  fields  are  described  using  Laplace  and  Poisson  equations 
given  as, 

v20  =  0,  V2Vr  =  -—  (8) 

££q 

where  £  is  the  dielectric  permittivity  and  e0  is  the  permittivity 
of  vacuum.  The  ion  transport  in  the  flow  field  is  governed  by  the 
extended  Nernst-Planck  equation  written  as, 


V  •  [DjVCj  +  m2-ZjCjV(0  +  VO  -  Vct\  =  0 


du  _  dv  _  dijf  _  d<fi  _  dq  _ 
dr  =  dr  =  !)F  =  lF  =  !F=u 


(12b) 


(9) 


Fig.  2.  Computational  domain. 
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Boundaries  CD  and  FE: 

dq  30  30 

U  =  0,  v  =  0,  —i  =  =  0 

dz  dz  dz 

Boundary  DE: 

30  30  a  dq 

u  =  0,v  =  0,  ^-  =  0,  -£-  = - ,  -^=0 

dr  or  ££q  or 

Boundary  GH: 

P2  =  0,  0  =  0,  02  —  0,  q  =  c0 


(12c) 

(12d) 

(12e) 


As  shown  in  Fig.  2,  boundaries  AB,  BC,  FG  and  GH  are  used 
to  emulate  the  infinitely  large  reservoirs.  Because  boundaries  AB 
and  GH  are  far  away  from  the  capillary  inlet  and  exit,  therefore 
the  pressure,  potential  and  ion  concentration  values  are  equal  to 
those  of  bulk  fluid  inside  the  reservoirs.  At  the  BC  and  FG  bound¬ 
aries,  the  fluid  velocity,  ion  concentration  and  potential  gradients 
in  radial  direction  are  specified  to  be  zero  to  emulate  the  infinite 
large  reservoirs.  The  CD  and  EF  boundaries  are  the  reservoir  walls.  It 
is  assumed  the  reservoir  walls  are  electrically  insulated  and  no  ion 
deposition  occurs.  For  the  boundary  condition  at  the  charged  wall, 
as  shown  in  Eq.  (12e),  we  assume  that  inside  the  Stern  layer  the 
ions  are  rigidly  held  and  do  not  contribute  to  the  ionic  distribution. 
The  surface  charge  density  satisfies  the  electroneutrality  condition 
at  the  wall. 


0.6 
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Fig.  3.  Comparison  of  maximum  streaming  potentials  computed  from  the  present 
numerical  model  with  the  results  of  Daiguji  et  al.  [  11  ]  under  the  same  computational 
domain  and  numerical  parameters. 


3.  Numerical  model 

Eqs.  (6)-(9)  along  with  the  boundary  conditions  given  in 
Eq.  (12)  were  solved  using  the  commercial  code  Comsol  3.3 
(http://www.comsol.com/).  The  finite  element  calculations  were 
performed  using  quadratic  triangular  elements.  Since  the  accuracy 
of  the  numerical  solutions  strongly  depends  on  the  mesh  size,  a 
refined  mesh  is  necessary  in  the  region  near  the  surface  where 
the  dependent  variable  gradients  are  pronounced.  A  finer  mesh 
is  used  near  the  charged  surface  to  capture  the  subtle  changes  in 
the  ion  concentration  and  electrical  potential.  Solution  indepen¬ 
dence  on  the  mesh  size  was  carefully  studied  before  reporting 
the  final  results.  The  numerical  results  show  that  the  solutions 
become  mesh-independent  when  the  element  number  exceeds 
approximately  6000.  Hence,  the  results  presented  in  this  study 
were  all  obtained  using  a  mesh  containing  14,000-20,000  ele¬ 
ments. 

4.  Results  and  discussion 

KC1  aqueous  solutions  at  temperature  T  =  298  K  were  used  as  the 
working  fluid.  The  dielectric  constant  of  the  KC1  aqueous  solution 
is  80,  the  molecular  diffusivities  of  K+  and  Cl-  ion  are  1.96  x  10-9 
and  2.03  x  10_9m2  s-1,  respectively.  The  density  and  viscosity  of 
the  solution  are  assumed  to  be  1000 kg m-3  and  10-3  Pas,  respec¬ 
tively.  The  ranges  of  bulk  concentration  c0,  capillary  radius  a  and 
surface  charge  density  a  are  chosen  as  10_5-10-2  M,  10-100 nm 
and  -0.25  x  10-4  to  -5  x  10-3  C  m-2,  respectively.  The  length  of  the 
capillary  is  fixed  with  L  =  5  [xm.  The  size  of  computational  domain 
in  the  reservoirs  is  fixed  as  Ar  x  Lr  =  1  |xm  x  1  |xm.  For  the  chosen 
capillary  radius  and  bulk  concentration  ranges,  the  corresponding 
dimensionless  Debye  length  Ka{Ka  =  a{££0RTl(2c0F2yi2)  is  ranged 
from  0.1  to  30. 

To  verify  the  correctness  of  the  proposed  numerical  model,  the 
maximum  streaming  potential  denoted  as  A0max  computed  from 
the  present  numerical  model  is  compared  with  that  reported  in  the 
study  of  Daiguji  et  al.  [11  ]  using  the  same  planar  slit  with  a  height 
of  30  nm.  The  boundary  conditions  for  the  planar  slit  are  the  same 
as  those  shown  in  Eq.  (12)  except  that  the  coordinate  system  is 
changed  from  (r,z)  to  (x,y).  As  shown  in  Fig.  3,  the  agreement  is  very 


good.  Based  on  this  comparison,  the  numerical  model  developed 
in  this  study  is  correct  and  can  be  extended  to  study  the  energy 
conversion  with  other  physical  domains.  In  Fig.  3,  A0max  for  the 
cylindrical  capillary  with  a  =  15  nm  is  also  computed  under  the  same 
parameters.  It  is  seen  that  the  variation  trend  of  A0max  is  similar  to 
that  of  planar  slit  but  with  smaller  magnitude.  This  implies  that  the 
energy  conversion  efficiency  predicted  using  the  cylindrical  cap¬ 
illary  would  be  smaller  than  that  using  the  planar  slit  since  the 
efficiency  is  proportional  to  A0max. 

As  discussed  in  the  one-dimensional  model,  the  maximum  effi¬ 
ciency  of  the  EI<  cell  is  independent  from  the  applied  pressure 
difference,  the  applied  pressure  difference  is  chosen  as  500  kPa 
in  this  study.  The  typical  results  of  ion  concentration,  electric 
potential,  velocity  field  and  current  flow  in  the  EK  cell  are  shown 
in  Figs.  4-7  for  various  streaming  potentials.  In  these  results, 
surface  charge  density  and  KCl  bulk  concentration  are  fixed  as 
o  =  -3  x  10-3  C  m-2  and  Co  =  10~4  M,  respectively.  As  shown  in  Fig.  4, 
the  ion  concentration  is  governed  mainly  by  the  surface  charge 
density  and  the  fluid  flow.  Large  variations  in  ion  concentration 
distribution  can  be  found  at  the  capillary  entrance  and  exit  regions. 
Inside  the  capillary,  concentrations  of  both  I<+  and  Cl-  remain  con¬ 
stants.  The  K+  ion  concentration  is  high  in  order  to  balance  the 
negative  charges  on  the  wall,  while  the  Cl-  ions  are  expelled  out 
of  the  capillary  and  have  a  much  lower  concentration  compared 
with  the  K+  ions.  Fig.  5  shows  the  combined  electrical  potential 
(0  +  0)  along  the  centerline  and  across  the  capillary  cross-section 
at  z  =  3.5  [xm.  As  shown  in  Fig.  5(a),  the  combined  electric  potential 
linearly  increases  from  the  inlet  to  outlet  reservoirs  and  the  fluid 
flow  does  not  affect  the  potential  distribution  significantly  except 
at  the  capillary  entrance  and  exit  because  of  geometric  change. 
For  the  A0  =  O  case,  electric  potential  can  still  be  found  because  of 
the  variation  in  electrostatic  potential  0.  Fig.  5(b)  shows  the  com¬ 
bined  electric  potential  profiles  across  the  capillary  cross-section  at 
z  =  3.5  [xm.  Because  the  potential  in  the  outlet  reservoir  is  chosen  as 
the  reference  and  the  surface  charge  density  is  negative,  the  com¬ 
bined  electric  potential  has  a  negative  value.  Note  that  the  potential 
profile  for  the  A0  =  O  case  is  due  to  electrostatic  potential  only. 
The  axial  velocity  along  the  capillary  centerline  shown  in  Fig.  6(a) 
becomes  a  constant  for  all  the  streaming  potentials  studied  indi¬ 
cating  the  flow  becomes  fully  developed  inside  the  capillary.  The 
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fully-developed  velocity  profiles  taken  at  z  =  3.5  p,m  are  shown  in 
Fig.  6(b).  Since  the  streaming  potential  results  in  an  electroosmotic 
flow  with  direction  opposite  to  the  pressure  driven  flow,  the  veloc¬ 
ity  is  reduced  as  the  streaming  potential  increases.  This  is  known  as 
the  electroviscous  effect.  In  Fig.  7(a),  the  current  density  along  the 
capillary  centerline  has  a  constant  value  inside  the  capillary  show¬ 
ing  that  the  current  flow  has  also  become  fully  developed.  The  fully 
developed  current  density  profile  at  z  =  3.5  p,m  shown  in  Fig.  7(b) 
indicates  that  the  current  flow  has  a  direction  opposite  to  the  fluid 
flow  in  the  region  near  the  wall.  In  this  region,  the  K+  ions  are  dom¬ 
inated  due  to  the  negative  surface  charge  and  the  fluid  velocity  is 
slow  due  to  the  no-slip  boundary  condition  at  the  wall.  Because  of 
the  electric  potential  distribution  shown  in  Fig.  5,  conductive  cur¬ 
rent  with  direction  opposite  to  the  streaming  current  is  larger  than 
the  streaming  current  and  results  in  negative  current  flow. 

Using  the  results  shown  in  Figs.  6  and  7,  one  can  construct  the 
electric  current- voltage  curve  ( I-V  curve)  and  flow  rate-voltage 
curve  (Q-  V  curve)  of  the  EI<  cell.  Fig.  8  shows  the  I-V  and  Q-  V  curves 


r/a 

Fig.  4.  Ion  concentration  distributions  in  EK  cell  under  various  streaming  potentials. 
Co  =  10-4  M,  a  =  50  nm,  cr  =  -3  x  10-3  C  rrr2,  A p  =  -500  kPa.  (a)  Along  the  centerline, 
(b)  Cross-section  at  z  =  3.5  (Jim. 
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Fig.  5.  Combined  electric  potential  (0-M/O  distributions  in  EI<  cell  under  various 
streaming  potentials.  Co  =  10-4  M,  a  =  50  nm,  cr  =  -3  x  10-3  C  m-2,  A p  =  -500  kPa.  (a) 
Along  the  centerline,  (b)  Cross-section  atz  =  3.5  |jim. 


for  various  bulk  concentrations  with  a  =  50  nm,  cr  =  - 3  x  10-3  C  m-2 
and  A p  =  -500  kPa.  The  electrical  current  and  flow  rate  are  obtained 
by  integrating  the  velocity  and  current  density  profiles  over  the 
capillary  cross-section  shown  in  Figs.  6  and  7(b).  From  the  results 
shown  in  Fig.  8,  the  I-V  and  Q-V  curves  are  linear  when  Co  =  10-3  M. 
As  c0  decreases  further  to  10-5  M,  nonlinear  I-V  and  Q-V  curves  are 
resulted.  As  shown  in  Fig.  8(a),  the  current  magnitude  is  relatively 
small  compared  with  other  bulk  concentrations  shown  in  Fig.  8(a). 
The  reason  for  this  result  will  be  discuss  later.  Similar  to  the  I-V 
curves,  nonlinear  Q-V  curve  is  found  when  Co  is  low  as  shown  in 
Fig.  8(b).  Fig.  8(b)  also  shows  that  the  effect  of  streaming  potential 
and  bulk  concentration  on  flow  rate  is  not  significant  since  the  flow 
rate  varies  in  a  small  range. 

From  the  I-V  and  Q-V  curves,  maximum  streaming  potential 
A0max»  maximum  current  /max  and  maximum  flow  rate  Qmax  can 
be  found.  In  Fig.  9,  A0max,  /max  and  Qmax  as  functions  of  c0  and  a 
under  the  conditions  of  A p  =  -500  kPa  and  a  =  50  nm  are  presented. 
Under  the  zero  current  condition,  A0max  can  be  found  and  shown  in 
Fig.  9(a).  It  is  seen  that  A0max  increases  with  the  increase  in  a  when 
c0  are  10-2  and  10-3  M.  When  c0  reduces  to  10-4  and  10-5  M,  A0max 
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is  found  to  increase  to  a  maximum  value  and  then  to  decrease  with 
the  increase  in  cr.  The  variation  trend  of  A0max  shown  in  Fig.  9(a) 
agrees  well  with  that  predicted  in  the  study  of  Daiguji  et  al.  [11]. 
/max  and  Qmax  are  obtained  under  the  A 0  =  0  condition.  Similar  to 
the  A0max,  /max  is  found  increase  with  the  increase  of  cr  when  c0  are 
equal  to  10-2  and  10-3  M  as  shown  in  Fig.  9(b).  When  c0  are  10-4 
and  10-5  M,  /max  is  found  to  increase  to  a  maximum  value  and  then 
to  decrease  with  the  increase  in  cr.  The  variation  in  Qmax  as  functions 
of  c0  and  cr  are  shown  in  Fig.  9(c).  For  cases  of  c0  equal  to  10-2  and 
10-3  M,  Qmax  remains  approximately  constant  as  a  increases.  For 
c0  equal  to  10-4  and  10-5  M,  Qmax  is  found  to  decrease  with  the 
increase  in  cr. 

The  physical  reason  for  the  results  shown  in  Fig.  9  was  explained 
in  detail  in  the  studies  of  Daiguji  et  al.  [11-12]  from  the  variation 
in  ionic  distribution  characteristics.  An  alternative  explanation  is 
given  from  the  variations  in  electric  potential  and  streaming  cur¬ 
rent  density  distributions  in  this  study.  As  shown  in  Fig.  10(a),  the 
electrical  potential  variation  along  the  capillary  exists  when  A 0  =  0. 


r/a 


Fig.  6.  Fluid  flow  velocity  variations  in  EK  cell  under  various  streaming  potentials. 
Co  =  10-4  M,  a  =  50  nm,  cr  =  -3  x  10-3  C  irr2,  A p  =  -500  kPa.  (a)  Along  the  centerline, 
(b)  Cross-section  at  z  =  3.5  (Jim. 


Fig.  7.  Electric  current  variations  in  EK  cell  under  various  streaming  potentials. 
Co  =  10-4  M,  a  =  50  nm,  a  =  -3  x  10-3  C  rrr2,  A p  =  -500  kPa.  (a)  Along  the  centerline, 
(b)  Cross-section  atz  =  3.5  p,m. 

Since  A0  is  zero,  this  potential  variation  is  obviously  contributed 
by  the  electrostatic  potential  which  is  produced  by  the  ion  con¬ 
centration  change  at  the  capillary  entrance  and  exit  and  the  fluid 
flow.  Fig.  10(a)  also  shows  that  the  electrostatic  potential  magni¬ 
tudes  and  its  gradient  increase  with  the  decrease  in  c0.  The  effect 
of  this  electrostatic  potential  gradient  is  to  produce  a  conductive 
electric  current  opposite  to  the  streaming  current.  Fig.  10(b)  shows 
the  electrical  current  density  profiles  across  the  capillary  cross- 
section  under  the  A</>  =  0  condition.  Because  of  the  electrostatic 
potential  gradient,  conductive  current  with  direction  opposite  to 
the  streaming  current  appears.  In  the  near-wall  region,  negative 
electric  current  is  found  because  of  strong  conductive  current  and 
weak  streaming  current. 

Note  that  the  electrostatic  potential  gradient  exists  not  only 
when  A0  =  O  but  also  in  the  non-zero  A0  cases.  This  implies  that 
the  existence  of  electrostatic  potential  gradient  in  the  flow  field 
will  reduce  the  current  magnitude,  especially  when  the  bulk  con¬ 
centration  is  low.  The  electrostatic  potential  gradient  appearance 
explains  the  significant  reductions  in  current  in  the  I-V curve  shown 
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in  Fig.  8(a)  and  maximum  current  /max  shown  in  Fig.  9(b)  when 
the  bulk  concentration  is  low.  From  the  streaming  potential  theory, 
maximum  streaming  potential  is  obtained  when  the  total  current 
is  zero.  Based  on  the  above  discussion  on  the  electric  current  char¬ 
acteristics,  nonlinear  variation  of  maximum  streaming  potential 
as  function  of  surface  charge  density  for  low  bulk  concentrations 
is  expected,  as  shown  in  Fig.  9(a).  In  the  one-dimensional  model, 
the  potential  gradient  due  to  the  electrostatic  potential  variation 
does  not  appear  since  the  capillary  is  assumed  to  be  infinitely 
long  and  without  considering  the  inlet  and  outlet  reservoirs.  Over¬ 
estimated  energy  conversion  is  expected  from  the  one-dimensional 
model. 

From  the  I-V and  Q-  Vcurves  discussed  above,  the  energy  conver¬ 
sion  efficiency  of  the  EI<  cell  can  be  computed  using  Eq.  (3 ).  In  Fig.  11 , 
the  EI<  cell  efficiencies  as  a  function  of  the  streaming  potential  under 
various  bulk  concentrations  are  shown.  When  A0/A0max  =  O  and 
A0/ A0max  =  1,  the  efficiencies  are  zero  because  of  zero  stream- 
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Fig.  9.  Maximum  streaming  potential  A0max,  maximum  electric  current  Jmax  and 
maximum  flow  rate  in  EI<  cell  as  functions  of  surface  charge  density  and  KC1  bulk 
concentration,  a  =  50  nm,  Ap  =  -500  kPa.  (a)  Streaming  potential,  (b)  Electric  current, 
(c)  Flow  rate. 


Fig.  8.  Electric  current-voltage  (I-V)  curves  and  flow  rate-voltage  (Q-V)  curves 
of  EI<  cell  for  various  KC1  bulk  concentrations.  a  =  50nm,  cr  =  -3  x  10-3  Cm-2, 
Ap  =  -500  kPa.  (a)  I-V  curves,  (b)  Q-V  curves. 
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ing  potential  and  zero  current,  respectively.  The  efficiency  profiles 
are  all  in  parabolic  shape  and  maximum  efficiency  depends  on  the 
bulk  concentration  and  streaming  potential.  For  the  case  of  c0  =  10-3 
and  10-4  M,  maximum  efficiency  occurs  at  about  A0/ A0max  =  0.5. 
For  c0  =  10-5  M,  the  maximum  efficiency  occurs  approximately  at 
A</>/  A0max  =  0.7  because  of  the  nonlinear  I-V  curve. 

The  results  shown  in  Fig.  11  indicate  that  the  maximum 
energy  conversion  efficiency  occurs  when  Co  =  10_4M,  a  =  50nm, 
a  =  -3  x  10-3  Cm-2,  Ap  =  -500kPa  and  A0/A</>max  =  0.5.  A  more 
detail  of  the  maximum  energy  conversion  efficiency  is  given  in 
Fig.  12  in  which  the  maximum  energy  conversion  efficiency  as 
functions  of  the  bulk  concentration,  surface  charge  density  and  cap¬ 
illary  radius  are  presented  and  compared  with  the  one-dimensional 
model  predictions  described  in  Eq.  (4).  The  dimensionless  Debye 
length  Ka  is  also  indicated  in  Fig.  12.  In  Fig.  12(a),  the  maximum 
efficiencies  with  c0  =  10-2  M  are  shown  for  capillary  radii  a  =  10,  50 
and  100  nm.  The  maximum  efficiency  increases  with  the  increasing 
surface  charge  density  and  decreasing  capillary  radius.  The  agree¬ 
ment  between  the  one-  and  two-dimensional  models  is  good  for  the 
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Fig.  10.  Combined  electric  potential  distribution  along  the  capillary  centerline  and 
current  density  profiles  at  z  =  3.5  |jim  as  function  of  KC1  bulk  concentration  under  the 
conditions  of  A0  =  0,  a  =  50  nm,  a  =  -3  x  10-3  C  irr2  and  A p  =  -500  kPa.  (a)  Potential, 
(b)  Current  density. 


Fig.  11.  Thermodynamic  efficiency  of  EK  cell  as  function  of  KC1  bulk  concentration 
and  streaming  potential,  a  =  50  nm,  a  =  -3  x  10-3  C  rrr2  and  A p  =  -500  kPa. 


surface  charge  density  and  capillary  radius  ranges  studied.  Small 
discrepancy  exists  for  the  case  of  a  =  10  nm.  For  the  c0  =  10-3  M  case 
shown  in  Fig.  12(b),  maximum  efficiency  also  increases  with  the 
increasing  surface  charge  density  and  decreasing  capillary  radius. 
Good  agreement  between  the  one-  and  two-dimensional  predic¬ 
tions  can  still  be  obtained  for  the  cases  of  a  =  50  and  100  nm 
(/ca  =  5.21  and  10.04)  for  the  surface  charge  density  range  stud¬ 
ied.  For  the  case  of  a  =  10nm  (/ca  =  1.04),  small  deviation  between 
the  one-  and  two-dimensional  models  appears.  For  the  case  of 
Co  =  10-4  M  shown  in  Fig.  12(c),  the  maximum  efficiency  increases 
to  a  maximum  value  and  then  decreases  slightly  as  the  surface 
charge  density  increases.  It  is  also  found  that  the  efficiency  for 
a  =  10nm  becomes  smaller  than  the  a  =  50  and  100  nm  cases.  The 
agreement  between  one-  and  two-dimensional  models  can  only 
be  found  when  /ca  is  large  and  surface  charge  density  is  small. 
For  small  /ca  and  large  surface  charge  density  cases,  deviation 
between  the  one-  and  two-dimensional  models  becomes  signifi¬ 
cant.  When  the  bulk  concentration  is  further  reduced  to  10-5  M  as 
shown  in  Fig.  12(d),  the  variation  trends  for  the  maximum  efficiency 
are  approximately  similar  to  those  found  for  c0  =  10-4  M  case  but 
with  smaller  magnitudes  and  more  significant  drops  in  the  high 
surface  charge  density  range.  The  predictions  between  one-  and 
two-dimensional  models  completely  disagree.  Since  the  assump¬ 
tions  made  in  the  one-dimensional  model  have  been  removed 
in  the  two-dimensional  model,  the  results  obtained  from  the 
two-dimensional  model  are  considered  to  be  more  realistic  in  pre¬ 
dicting  the  electrokinetic  energy  conversion  performance.  From 
the  results  shown  in  Fig.  12,  it  is  seen  that  the  one-dimensional 
model  is  valid  only  for  large  /ca  and  small  surface  charge  den¬ 
sity. 

Since  the  efficiency  is  evaluated  based  on  the  variations  in  cur¬ 
rent,  streaming  potential  and  flow  rate,  we  use  the  results  for  the 
a  =  50  nm  case  to  explain  the  observations  shown  in  Fig.  12  in  more 
detail.  In  Fig.  9,  typical  results  of  maximum  streaming  potential, 
current  and  flow  rate  are  shown.  Since  the  applied  pressure  differ¬ 
ence  is  fixed  as  A p  =  -500  kPa  and  flow  rate  remains  approximately 
constant,  the  maximum  efficiency  would  be  proportional  to  the 
magnitude  of  maximum  power  output  A</>max/max.  For  c0  =  10-2 
and  10-3  M  cases,  both  A0max  and  /max  increase  with  the  surface 
charge  density.  Therefore,  the  maximum  efficiency  increases  with 
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Fig.  12.  Maximum  efficiencies  as  functions  of  KC1  bulk  concentration,  capillary  radius  and  surface  charge  density,  (a)  Co  =  10-2  M.  (b)  Co  =  10-3  M.  (c)  Co  =  10-4  M.  (d)  Cq  =  10-4  M. 


the  increase  in  surface  charge  density  as  indicated  in  Figs.  12(a  and 
b).  For  Co  =  10-4  and  10-5  M  cases,  variations  of  A0max  and  /max 
become  nonlinear  with  respect  to  the  surface  charge  density.  The 
resulting  maximum  efficiency  also  varies  nonlinearly  with  the  sur¬ 
face  charge  density.  Because  of  the  decreases  in  A0max  and  Jmax 
when  the  surface  charge  density  is  large,  the  maximum  efficiencies 
for  these  two  bulk  concentrations  also  decrease  when  the  surface 
charge  is  high,  as  shown  in  Fig.  12(c  and  d). 

Fig.  13  shows  the  maximum  efficiency  as  a  function  of  /ca  for 
various  surface  charge  densities  to  provide  better  understanding 
of  the  maximum  efficiency  dependence  on  /ca.  It  is  seen  that  an 
optimum  maximum  efficiency  occurs  when  /ca  is  equal  to  2  and  its 
magnitude  increases  with  the  increase  in  surface  charge  density. 

Besides  the  EK  cell,  another  important  application  of  electroki- 
netic  effect  is  the  electrokinetic  pump  (EK  pump)  [20,21  ].  Since  the 
EK  pump  operation  is  the  reverse  of  the  EK  cell,  the  current  numer¬ 
ical  model  can  be  easily  extended  to  the  2D  simulation  of  EK  pump. 
In  EK  pump,  the  fluid  flow  is  driven  by  an  applied  voltage  and  a 


pressure  rise  is  created  in  the  flow  direction.  Similar  to  the  EK  cell 
case,  EK  pump  operation  can  be  described  by  using  the  pressure  and 
voltage  differences  between  the  inlet  and  outlet  reservoirs.  Refer¬ 
ring  to  Fig.  1 ,  pressure  and  potential  in  the  inlet  and  outlet  reservoirs 
are  denoted  as  0i ,  Pi ,  02,  P2>  respectively.  Using  the  current  estab¬ 
lished  model  to  simulate  the  EK  pump,  an  externally  applied  voltage 
A0  =  0 2  -0i  <0  and  a  back  pressure  Ap=p2-pi  >0  are  specified 
at  the  inlet  and  outlet  reservoirs.  Referring  to  Fig.  2,  the  boundary 
conditions  needed  to  be  changed  are 

Boundary  AB: 

Pi  =  0,  0  =  0,  0i  =  -  A 0,  C/  =  c0 
Boundary  GH: 

p2  =  A p,  02  =  0,  0  =  0,  q  =  c0 


470 


R.  Chein  et  al.  /  Journal  of  Power  Sources  187  (2009)  461-470 


Fig.  13.  Maximum  efficiency  as  functions  of  dimensionless  Debye  length  Ka  for 
various  surface  charge  densities. 


5.  Conclusion 

The  energy  conversion  efficiency  of  electrical  power  genera¬ 
tion  using  the  electrokinetic  effect  was  investigated  numerically.  A 
two-dimensional  axisymmetric  physical  domain  including  a  finite- 
length  charged  cylindrical  capillary  and  reservoirs  connected  at 
the  capillary  ends  was  considered.  The  fluid  flow  velocity,  ion 
concentration  and  electrical  potential  were  solved  simultaneously 
based  on  Navier-Stokes,  Nernst-Planck,  Laplace  and  Poisson  equa¬ 
tions.  Using  the  computed  velocity,  ion  concentration  and  electrical 
potential  distributions,  the  electrical  current  flow  in  the  capillary 
was  computed. 

For  the  bulk  concentration,  capillary  radius  and  surface  charge 
density  in  the  ranges  of  10_5-10-2  M,  10-100  nm  and  -0.25  x  10~4 
to  -5  x  10-3  Cm-2,  the  two-dimensional  computational  results 
show  that  the  electrostatic  potential  gradient  produced  by  the 
changes  in  ion  concentration  at  the  capillary  entrance  and  exit  and 
fluid  flow  plays  an  important  role  in  evaluating  EK  cell  efficiency. 
It  was  found  that  the  electrostatic  potential  gradient  produces  a 
conductive  current  in  direction  opposite  that  to  the  streaming 
potential,  in  addition  to  that  produced  by  the  streaming  poten¬ 
tial.  It  was  also  found  that  the  conductive  current  produced  by 
the  electrostatic  potential  gradient  increased  with  the  decrease 
in  bulk  concentration  and  increase  in  surface  charge  density.  In 
terms  of  the  dimensionless  Debye  length  /ca,  this  indicated  that 
the  electrostatic  potential  gradient  is  large  when  /ca  is  small.  In  the 


one-dimensional  model,  the  electrostatic  potential  gradient  can¬ 
not  be  predicted  because  the  capillary  inlet  and  exit  are  excluded. 
This  implies  that  the  electrical  current  predicted  from  the  one¬ 
dimensional  model  would  be  higher  than  that  predicted  from  the 
two-dimensional  model. 

Because  of  the  conductive  current  produced  by  the  electrostatic 
potential  gradient,  the  total  current  is  reduced  in  the  flow  field. 
The  computed  results  show  that  current-voltage  curve  ( l-V  curve) 
and  flow  rate-voltage  curve  (Q-V  curve)  become  nonlinear  when 
the  bulk  concentration  is  low.  Using  the  l-V  and  Q-V  curves,  the 
maximum  EK  cell  efficiency  is  evaluated.  For  high  bulk  concentra¬ 
tion  cases,  it  was  found  that  maximum  efficiency  increased  with 
the  increase  in  surface  charge  density.  For  low  bulk  concentra¬ 
tion  cases,  the  maximum  efficiency  increased  to  a  maximum  value 
and  then  decreased  slightly  with  respect  to  the  increase  in  surface 
charge  density.  The  comparison  between  the  maximum  efficien¬ 
cies  predicted  from  one-  and  two-dimensional  models  indicated 
that  the  one-dimensional  model  is  valid  only  when  /ca  is  large  and 
the  surface  charge  density  is  low.  For  small  /ca  and  large  surface 
charge  cases,  the  one-dimensional  model  over-predicted  the  effi¬ 
ciency  because  the  electrostatic  potential  gradient  is  neglected. 
An  optimum  maximum  efficiency  is  found  when  /ca  is  equal  to  2 
and  its  magnitude  increases  with  the  increase  in  surface  charge 
density  based  on  the  two-dimensional  model  results. 
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